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We present a model that elucidates the physics underlying 
the generation of an electromagnetic pulse from a femtosec- 
ond laser induced plasma channel. The radiation pressure 
force from the laser pulse spatially separates the ionized elec- 
trons from the heavier ions and the induced dipole moment 
subsequently oscillates at the plasma frequency and radiates 
an electromagnetic pulse. 
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Recent investigations of the propagation of intense 
femtosecond infrared (IR) laser pulses in air show that 
the dynamical interaction between nonlinear self-focusing, 
plasma defocusing, and group-velocity dispersion can 
cause an initial beam to break up spatially into sev- 
eral filaments, or light strings, with diameters around 
a hundred microns that can maintain themselves over 
long distances [|]-|| . It has been observed experimentally 
that femtosecond light strings in turn produce plasma 
channels by multi-photon ionization (MPI) along their 
direction of propagation with lengths on the order of 
centimeters ]6|-|i"(| . Observations of the electromagnetic 
pulses (EMPs) from light string induced plasmas suggest 
that these channels attain dipole moments during the 
laser pulse which subsequently oscillate at the plasma 
frequency and radiate ]l0|| . 

In this letter we present a simple model for EMP gen- 
eration from femtosecond light string induced plasmas to 
elucidate the underlying physics, in particular how the 
initial dipole moment is produced in the plasma chan- 
nel, and the characteristic field strengths and spectrum 
of the EMP. In this model electrons and ions are initially 
created by MPI due to the light string, and a dipole mo- 
ment is subsequently induced in the plasma via the radia- 
tion pressure force resulting from the velocity-dependent 
Lorentz force. For a laser field of frequency u> this force 
contains a second-harmonic component at 2uj in addi- 
tion to a low frequency component, and it is this low 
frequency component of the force that creates the initial 
dipole moment in the plasma channel by longitudinally 
separating the light electrons from the heavy ions on the 
short time scale of the light string. Following the laser 
pulse the electron-ion system then starts to oscillate lon- 
gitudinally at the plasma frequency and generates an 
EMP propagating radially away from the plasma chan- 
nel Jl0[] . We stress that this mechanism for creating a 
dipole moment in a plasma is new and distinct from that 



due to the ponderomotive force from a tightly focused 
Gaussian beam 12 r|: In particular, the usual pondero- 
motive force acting on electrons is proportional to the 
spatial gradient of the applied light intensity jl4j]. For 
light strings of diameter 100 /im and a few centimeters 
length the ponderomotive force should be negligible, and 
this was our motivation for looking for a new physical 
process that could create a dipole moment in the plasma 
channel. 

In our model we consider a femtosecond laser pulse 
linearly polarized along the x-axis, and propagating in air 
along the positive z-axis, with a central frequency lu in the 
IR region. For an electron plasma in the presence of an 
electric field E and a magnetic field B, the hydrodynamic 
equation for the average velocity of an electron v reads 
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where e and m are the charge and mass of the electron, 
and (dv/dt)\ scaU is the time rate of change of the velocity 
due to constituent particle scattering. The parameter rj 
is included to keep track of orders of perturbation theory 
and reflects the fact that the velocity-dependent Lorentz 
force is order \v\/c << 1 smaller than the dominant force 
due to the electric field under the non-relativistic condi- 
tions appropriate here. To keep the theory as simple as 
possible we concentrate on the electron motion, assum- 
ing that the heavy ions remain stationary. In addition, 
we have the continuity equation for the electron density 
I5H 
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where N e is the electron density, u> is the central fre- 
quency of the input laser pulse, K is the order of the 
MPI and the sum allows for the various constituents of 
air, mainly oxygen and nitrogen molecules, (3^ K ^ is the 
MPI coefficient derived from, for example, Keldysh the- 
ory jlffl , and I{r,t) is the intensity distribution of the 
laser pulse. The right-hand side of Eq. (||) describes 
the generation of electrons via MPI, since the effects of 
avalanche ionization for femtosecond pulses are negligi- 
ble compared to MPI. To complete the system we couple 
these equations to the Maxwell equations for the elec- 
tromagnetic field incorporating the effects of nonlinear 
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self-focusing, plasma defocusing, and group velocity dis- 
persion. The key Maxwell equation we need here is 



V • E = -e{N e - Ni)/e , 



(3) 



up ~ 10 fs. Due to temporal compression associated with 

, the duration t p of the 



which determines the irrotational part of the electric field 
from the net charge density —e(N e — Ni)/eo, with Ni the 
ion density of the plasma. 

Equations provide the basic dynamics of the 

plasma. Detailed numerical simulations of the formation 
and structure of light strings have been performed in a 
variety of situations Typically, light strings have 

a diameter of around 100 /im and a temporal duration 
t 

sclf-focusing of the input pulse 
light string is much shorter than the input pulse duration. 
Here we shall use electric and magnetic fields which re- 
flect the characteristic properties of a single light string 
rather than explicitly solving for their detailed structure, 
which would obscure the physics to be presented. 

To proceed we first consider the evolution of the elec- 
tron velocity for short times on the order of the light 
string duration t p during which time the electrons are 
generated. For such short times we may safely neglect 
the effects of particle diffusion, hence we ignore spatial 
derivatives in Eq. (|l|). When the electrons are liberated 
via MPI they are accelerated in the laser field and un- 
dergo collisions, and we approximate the scattering as 
(dv/dt)\ scatt w -7«, with 7 « 10 12 - 10 13 s" 1 the colli- 
sion rate for air at atmospheric pressure ]lO||lS| ] . Further- 
more, the dominant electromagnetic fields will be those 
due to the light string E — El and B — Bl- The spa- 
tial dimensions of the light string and plasma channel are 
around 100 /im 2 x 1 cm, enormous on the scale of a cubic 
wavelength, and we approximate the fields as spatially 
uniform over the transverse profile of the plasma chan- 
nel. The fields El and Bl due to the light string can 
then be expressed as 



E L {r,t) = -[£ L {t-z/v g y 



+ c.c 



(4) 



with B L {r,t) = yE L /c, where £ l (t) = y/2I /ce Q -e T ' ft p 
is the Gaussian approximation to the light string electric 
field envelope with t = t — z/v g the retarded time in a 
reference frame moving at the light string group velocity 
v g , Iq and t p being the peak intensity and duration of 
the light string, no ~ 1 is the refractive-index of air, 
and k = uj/c is the magnitude of the optical wave vector. 
Since the light string is polarized along the x-direction we 
expect the electron motion to be predominantly along the 
x-axis, and we correspondingly write the electron velocity 
as v = xv x + r\'zv z : By symmetry the electron motion will 
be confined to the x-z plane. Then to zeroth and first 
orders in the perturbation parameter i] Eq. (|l|) yields 
the pair of equations 
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Under the assumption that the electric field envelope 
£l{t) varies slowly on the carrier time scale w _1 the 
equation for v x may be solved approximately as 

v x (z,t)^~— M 1 9) + c.c . 6 

2m \ (7 — iuj) J 

Next, inspection of the driving force term proportional 
to ~e(v x j 'c)El appearing in Eq. (|^) for v z reveals that 
it has high frequency or second-harmonic terms varying 
as ±2w in addition to a low frequency term. The low 
frequency term produces a force that is always directed 
along the positive z-axis, and creates a charge separa- 
tion between the fast responding electrons and slow ions 
along the direction of the light string, and hence an ini- 
tial dipole moment. In contrast, the charge separation 
induced by the second-harmonic term keeps switching di- 
rection, hence producing zero dipole moment on average, 
and we therefore isolate the low frequency component. 
Then by substituting Eqs. (Q) and ^) into (||) for v z we 
obtain 
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where v z = d£/dt, £,(z,t) being the longitudinal (z) 
electron-ion separation, and Frp is the radiation pres- 
sure force exerted on the electrons (see below). Inte- 
grating Eq. (Q) over time and using £(z, —00) = and 
d£,/dt\±oo — we obtain 
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where / = ^qc\£l\ 2 is the light string intensity, we have 
used uj >> 7 to simplify the expression, and we have 
extended the upper time limit to infinity with the under- 
standing that we formally only need a large enough time 
window to cover the light string duration. 

The quantity £0 is the longitudinal electron-ion sepa- 
ration induced by the light string, with associated dipole 
moment directed along the z-axis do = — e£o, and is pro- 
portional to the time-integrated intensity or fluence of 
the light string, but relatively insensitive to the scatter- 
ing rate so long as uj >> 7 > 0. In general, £0 could 
depend on position along the light string due to longi- 
tudinal variation of the intensity profile, but in keeping 
with our approximations we shall treat £ as a constant 
over the length of the plasma channel. Physically, the 
low frequency component of the Lorentz force exerts a 
force on each electron in the direction of propagation of 
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the light string. To elucidate the physical origin of this 
force we realize that expression (^) for the radiation pres- 
sure force can be written as Frp — hk ■ (I/hui) ■ a, where 
hk is the momentum transferred to the electron along 
the z-axis per absorbed photon, (I/hui) is the photon 
flux (number of photons per unit area per unit time), 
and a — (fce 2 7/wme )/(7 2 + uj 2 ) is the cross-section 
for inverse Bremsstrahlung pj| . The force that causes 
the electron-ion separation along the direction of prop- 
agation of the light string may therefore be interpreted 
as the radiation pressure force that accompanies inverse 
Bremsstrahlung in which the electrons absorb photons in 
the field of the ions. 

Our physical picture of the EMP generation is now as 
follows: The light string of duration t p w 10 fs generates 
an initial dipole moment along the z-axis in the plasma 
channel which subsequently oscillates at the plasma fre- 
quency uipi — y/N e e 2 /me a Jll],^0| , and radiates an EMP 
radially outwards from the channel (see below). For elec- 
tron densities N e 10 15 cm~ 3 typical of light string 
induced plasmas, we have ui p i = 2 x 10 12 rads -1 with 
corresponding EMP wavelength Xemp = 27rc/ui p i ~ 1 
mm. Some other physical considerations are in order re- 
garding this picture. In particular, in addition to the 
longitudinal electron-ion displacement £o induced by the 
light string the electrons also attain temperatures of the 
order of a few electron-volts due to the oscillatory 
force —eE^ acting on them along the x-axis in the field 
of the light string. Physically, electron-electron collisions 
will cause the electron plasma, initially described by a 
non-equilibrium electron distribution, to relax towards 
an Maxwellian distribution which as a whole executes 
plasma oscillations along the z-axis |2l]]. We do not ex- 
pect plasma oscillations along the x-axis as the diameter 
of the plasma channel (~ 100 /im) is smaller than the 
EMP wavelength. The ionized electron-ion pairs can also 
undergo radiative recombination hence producing high 
frequency EM fields in addition to the relatively low fre- 
quency EMP of interest here. 

To calculate the plasma-electron oscillations of the 
channel for times after the light string has passed we 
use the classic approach of Tonks and Langmuir [p"i"|| : 
Then Eq. (||) for the self-consistent electric field re- 
sulting from the charge separation becomes dE z /dz — 
(eN e /e a )(d£/dz) 1 with solution E z = (eiV e /e )£. The 
longitudinal force on each electron is then —eE z — 
m<9£ 2 /dt 2 , which yields the following equation of motion 
for the electron-ion displacement 
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where oj p i{t) = \J N e {t)e 2 /meo as before. Since the 
plasma density varies on a length scale L of centimeters 
which is much larger than the anticipated EMP wave- 
length Xemp ~ 1 mm, we hereafter treat the plasma 



density, and hence the plasma frequency as constant over 
the plasma channel. Equation (^) is to be solved with 
the initial condition £(0) = £o from Eq. (|J): This as- 
sumes that the initial electron-ion separation is not sig- 
nificantly degraded by thermalizing collisions, meaning 
that we are taking an upper bound for £(0). For long 
times well after the light string has passed, the elec- 
tron density decays due to radiative recombination as 
dN e /dt = —aN 2 which has the analytic solution 
N e (t) = N e (0)/(1 + atN e {0)). The initial value of the 
electron density can be calculated by integrating Eq. (||) , 
but here we use typical value of N e (0) — 10 15 cm~ 3 ob- 
tained from the full numerical solutions 0-B. 

To obtain an approximate solution for the EMP due to 
the oscillating electron plasma we write the longitudinal 
current as J(r,t) = —eN e (d£/dt)z. Furthermore, since 
the plasma channel radius tq ~ 50 /im is smaller than the 
EMP wavelength A em p ~ 1 mm, and the plasma channel 
length L >> Xemp, we approximate the induced current 
in the channel as a line source 



J{r,t) » -eN e {t)i{t)8{x)5{y)^r 2 )z 1 
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where the dot signifies a time derivative. This current os- 
cillating along the z-axis will radiate a radially symmetric 
EMP E ra d(f,t) = zE ra d(r,t), independent of z and de- 
pendent only on the radial distance r, that propagates 
through air and radially away from the plasma channel, 
and obeys Maxwell wave equation 
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It is important in this model that the EMP radiates into 
the air surrounding the plasma channel as radiation at 
the plasma frequency u> p i has zero group velocity and can- 
not therefore propagate in the plasma channel. Equation 
( |TT| ) along with the current in Eq. ([h]) can be solved 
using the Green's function for the wave equation [E2j 



E ra d(r, t) 



en N e (t)rl 
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where we have taken the density N e (t) outside the inte- 
gral based on the fact that it decays on a time scale of 
nanoseconds, slow compared to the plasma oscillations. 
Equation (^|) shows that even though we have not in- 
cluded damping of the plasma oscillation the radiated 
EMP is limited in duration by the time it takes the elec- 
tron plasma to recombine, which is on the order of a few 
nanoseconds [pHLOl: In contrast, Landau damping p0| of 
the plasma oscillations for our parameters occurred on a 
much longer time scale. 

To illustrate the features of the EMP we have solved 




and (|12|) numerically with the initial condi- 
Here we chose the light string parameters 
2 x 10 13 W/cm 2 ,t p = 10 fs, an initial electron den- 
sity iV e (0) = 10 15 cm- 3 , and a = 5 x 10~ 7 cm 3 /s §. We 



3 



find that the radiated EMP has a peak field strength of 
l-EWioxl ~ 2 kV/m at an observation distance r — 1 mm 
from the channel. The spectrum of the EMP is shown in 
Fig. (|l|), and the peak of the spectrum is at uj pea k ~ 10 12 
rad s -1 , which lies in the far IR region, which is con- 
sistent with the plasma frequency for a density around 
N e w 10 15 cm" 3 . The Hertzian width of the EMP spec- 
trum is around Av 3x 10 9 s -1 , which is consistent with 
the inverse lifetime due to radiative electron recombina- 
tion aN e {0) » 10 9 s _1 . Thus, the peak EMP frequency 
serves as a measure of the plasma density, and the width 
of the spectrum is a measure of the decay time of the 
electron plasma. This suggests that the EMP from light 
string induced plasma could serve as a useful diagnostic. 
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FIG. 1. Spectrum of the EMP. 

In summary, we have shown that the EMP from a light 
string induced plasma channel can be understood as due 
to the establishment of a dipole moment in the channel 
by the radiation pressure force of the light string, which 
subsequently oscillates at the plasma frequency and ra- 
diates. Having established the basic physics, in future 
work we shall explore the details of the evolution of the 
plasma channel using the full kinetic theory ]23jj24| ], and 
in particular how the initial dipole moment dephases be- 
fore radiating. This raises the prospect that the EMP 
field strength may be enhanced by creating M parallel 
plasma channels || within the EMP wavelength, then 
the EMP intensity should be enhanced by M 2 . 
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